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Abstract 

Lattice coding techniques may be used to derive achievable rate regions which outperform known 
independent, identically distributed (i.i.d.) random codes in multi-source relay networks and in particular 
the two-way relay channel. Gains stem from the ability to decode the sum of codewords (or messages) 
using lattice codes at higher rates than possible with i.i.d. random codes. Here we develop a novel lattice 
coding scheme for the Two-way Two-relay Channel: 1 2 «-» 3 O 4, where Node 1 and 4 simultaneously 
communicate with each other through two relay nodes 2 and 3. Each node only communicates with its 
neighboring nodes. The key technical contribution is the lattice-based achievability strategy, where each 
relay is able to remove the noise while decoding the sum of several signals in a Block Markov strategy and 
then re-encode the signal into another lattice codeword using the so-called "Re-distribution Transform". 
This allows nodes further down the line to again decode sums of lattice codewords. This transform is 
central to improving the achievable rates, and ensures that the messages traveling in each of the two 
directions fully utilize the relay's power, even under asymmetric channel conditions. All decoders are 
lattice decoders and only a single nested lattice codebook pair is needed. The symmetric rate achieved 
by the proposed lattice coding scheme is within ~ log 3 bit/Hz/s of the symmetric rate capacity. 

I. Introduction 

Lattice codes may be viewed as linear codes in Euclidean space; the sum of two lattice codewords is 
still a codeword. This group property is exploited in several additive white Gaussian noise (AWGN) 
relay networks, such as (lJ-Q, where it has been shown that lattice codes may outperform i.i.d. 
random codes in certain scenarios, particularly when interested in decoding a linear combination of 
the received codewords rather than the individual codewords. One such example is the AWGN Two-way 
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Relay channel, where two users wish to communicate with each other through a relay node (3), j4). 
If the two users employ properly chosen lattice codewords, the relay node may decode the sum of the 
codewords from both users directly at higher rates than decoding them individually. It is then sufficient 
for the relay to broadcast this sum of codewords to both users since each user, knowing the sum and its 
own message, may determine the other desired message. 

Past work. Lattice codes have been used to derive achievable rates for multi-user networks beyond 
the two-way relay channel. Nested lattice codes have been shown to be capacity achieving in the point- 
to-point Gaussian channel |5}, the Gaussian Multiple-access Channel |2], Broadcast Channel [6], and 
to achieve the same rates as those achieved by i.i.d. Gaussian codes in the Decode-and-Forward rate 
Q), and Compress-and-Forward rates [1] of the Relay Channel |8). Lattice codes may further be 
used in achieving the capacity of Gaussian channels with interference or state known at the transmitter 
(but not receiver) j9j using a lattice equivalent [6| of dirty-paper coding (DPC) [10]. The nested lattice 
approach of [6| for the dirty-paper channel is extended to dirty -paper networks in [11]. Lattice codes 
have also been shown to approximately achieve the capacity of the K-user interference channel fT2) , 
using the Compute-and-Forward framework for decoding linear equations of codewords (or equivalently 
messages) of Q. 

We consider the Two-way Two-relay Channel: 1 f> 2 f> 3 f> 4 where two user Nodes 1 and 4 
exchange information with each other through the relay nodes 2 and 3. This is related to the work of 
fl3) , which considers the throughput of i.i.d. random code-based Amplify-and-Forward and Decode-and- 
Forward approaches for this channel model, or the i.i.d. random coding based schemes of fl4[ and fl5) 
where there are furthermore links between all nodes. This model is also different from that in [16| where 
a two-way relay channel with two parallel (rather than sequential as in this work) relays are considered. 

Contributions. The Two-way Two-relay Channel 1 f> 2 f> 3 f> 4 is a generalization of the Two- 
way Relay channel 1 <H- 2 <H> 3 to multiple relays. As lattice codes proved useful in the Two-way Relay 
channel in decoding the sum at the relay, it is natural to expect lattice codes to again perform well for the 
considered channel. We propose, for the first time, a lattice based scheme for this two-way line network, 
where all nodes transmit lattice codewords and each relay node decodes a sum of these codewords. This 
scheme may be seen as a generalization of the lattice based scheme of (3), Q for the Two-way Relay 
Channel. However, this generalization is not straightforward due to the presence of multiple relays and 
hence the need to repeatedly be able to decode the sum of codewords. One way to enable this is to have 
the relays employ lattice codewords as well - something not required in the Two-way Relay channel. In 
the Two-way Relay channel achievability schemes consists of two phases - the multiple access and the 
broadcast phase. The scheme includes multiple Block Markov phases, where during each phase, the end 
users send new messages encoded by lattice codewords and the relays decode a combination of lattice 
codewords. The relays then perform a "Re-distribution Transform" on the decoded lattice codeword 
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combinations, and broadcasts the resulting lattice codewords. The novelty of our scheme lies in this 
"Re-distribution Transform" which enables the messages traveling in both directions to fully utilize the 
relay transmit power. Furthermore, all decoders are lattice decoders (more computationally efficient than 
joint typicality decoders) and only a single nested lattice codebook pair is needed. 

Outline. We first outline some definitions, notation, and technical lemmas for nested lattice codes in 
Section [n] We outline the channel model in Section [III] We then outline a lattice based strategy for the 



broadcast phase of the Two-way (single) Relay channel with asymmetric uplinks in Section IV which 
includes the key technical novelty - the "Re-distribution Transform" which allows relays to intuitively 
spread the signals traveling in both directions to utilize the relay's entire transmit power. We present the 
main achievable rate regions for the Two-way Two-relay channel in Section [VJbefore outlining extensions 



to half-duplex nodes and more than two relays in Section VI Finally, we conclude in Section VII 



II. Preliminaries on lattice codes and notation 



We now define our notation for lattice codes in Subsection II-A define properties of nested lattice 



codes in Subsection III-BI and technical lemmas in Subsection IH-CI 



A. Lattice codes 

Our notation for (nested) lattice codes for transmission over AWGN channels follows that of |6|, 
p7) ; comprehensive treatments may be found in (5), (5), and in particular |19|. An ri-dimensional 
lattice A is a discrete subgroup of Euclidean space W 1 with Euclidean norm 1 1 • 1 1 under vector addition 
and may be expressed as all integral combinations of basis vectors bj £ R™ 

A = {A = B i : i e Z"}, 

for Z the set of integers, n G Z + , and B := [bi |b2 1 ■ ■ • b n ] the n x n generator matrix corresponding 
to the lattice A. We use bold x to denote column vectors, x T to denote the transpose of the vector 
x. All vectors lie in R" unless otherwise stated, and all logarithms are base 2. Let denote the all 
zeros vector of length n, I denote the n x n identity matrix, and A/(//, a 2 ) denote a Gaussian random 
variable (or vector) of mean \i and variance a 2 . Let \C\ denote the cardinality of the set C. Define 
C(x) :— \ log 2 (1 + x). Further define or note that 

• The nearest neighbor lattice quantizer of A as 

Q(x) = argmin||x - A||; 
AeA 

• The mod A operation as x mod A := x — Q(x); 

• The Voronoi region of A as the points closer to the origin than to any other lattice point 

V := {x : Q(x) = 0}, 
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which is of volume V := Vol(V) (also sometimes denoted by V(A) or Vi for lattice Aj); 

• The second moment per dimension of a uniform distribution over V as 

• For any s G R n , 

(a(s mod A)) mod A = (as) mod A, a G Z. (1) 
/3(s mod A) = (/3s) mod /3A, /3 G M. (2) 

• The definitions of Rogers good and Poltyrev good lattices are stated in pi; we will not need these 
definitions explicitly. Rather, we will use the results derived from lattices with these properties. 

B. Nested lattice codes 

Consider two lattices A and A c such that A C A c with fundamental regions V, V c of volumes V, V c 
respectively. Here A is termed the coarse lattice which is a sublattice of A c , the fine lattice, and hence 
V > V c . When transmitting over the AWGN channel, one may use the set Ca c .v = {A c n V(A)} as the 
codebook. The coding rate R of this nested (A, A c ) lattice pair is defined as 

R= ilog |C Ac ,v| = -l°g 7^- 
n n V c 

Nested lattice pairs were shown to be capacity achieving for the AWGN channel |5|. 

In this work, we only need one "good" nested lattice pair A C A c , in which A is both Rogers 
good and Poltyrev good and A c is Poltyrev good (see definitions of these in [ 1 1). This lattice pair is 
used throughout the paper. The existence of such a pair of nested lattices may be guaranteed by [5|; we 
now describe the construction procedure of such a good nested lattice pair. Suppose the coarse lattice 
A = BZ" is both Rogers good and Poltyrev good with second moment er 2 (A) = 1. With respect to this 
coarse lattice A, the fine lattice A c is generated by Construction A Q, Q, which maps a codebook of 
a linear block code over a finite field into real lattice points. The generation procedure is: 

• Consider the vector G G Fp with every element drawn from an i.i.d. uniform distribution from 
the finite field of order P pr ime (a prime number) which we take to be Fp jme = {0, 1, 2, ... , P pr ime— 
1} under addition and multiplication modulo P pr ime- 

• The codebook C of the linear block code induced by G is C — {c — Gw : w G Fp prime }. 

• Embed this codebook into the unit cube by scaling down by a factor of P pr ime and then place a 
copy at every integer vector: A c = P pr l me C + Z™. 

• Rotate A c by the generator matrix of the coarse lattice to obtain the desired fine lattice: A c = BA C . 
Now let (/)(■) denote the one-to-one mapping between one element in the one dimensional finite field 
w G Fp primi , to a point in n-dimension real space t G Ca c ,v : 

t = 4>{w) = (BP-] me Gw) mod A, (3) 
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with inverse mapping w = <f) = (G T G) 1 G T (P p i rme (B 1 t mod Z n )) (see Q Lemma 5 and 
6]). The mapping operation <j>{-) defined here is used in the following lemmas. 

C. Technical lemmas 

We first state several lemmas needed in the proposed two-way lattice based scheme. In the following, 

t ai and tu e C Ac ,v are generated from w ai and w bi € Wp prirne as t ai = <j>(w ai ),t bl = <j>(w bi ). 
Furthermore, let a, on, (3. t e Z such that p-^— , -p^ 1 —, -p^— i Z and G K. We use 0, © and 
G to denote modulo P pr ime addition, multiplication, and subtraction over the finite field Fp primc . 

Lemma 1: There exists an one-to-one mapping between v = Q]\ cxi&t a i+^2 i Pi&bbi) mod OA and 

Proof: The proof follows from [2 Lemma 6], where it is shown that there is an one-to-one mapping 
between _1 v = (Y /i a i t ai + Z)< A*w) mod A and u' = 4 a -w ai © 4 /3> fci , where a- = a 4 
mod Pprime and = mod P pr ime- This one-to-one mapping is given by _1 v = ^ _1 (u') and 
u' = 0(0 _1 v). Then observe that 0^ ajuw © 4 /3iW bi = i o' l w ai © l /3;U>bi by the properties 
of modulo addition and multiplication []. Thus, the one-to-one mapping between v and u is given by 
v = 00" 1 (u) and u = 0(0 -1 v). ■ 
Lemma 2: There exists an one-to-one mapping between a © u> and u>. 

Proof: Recall that -p-^ — ^ Z. Suppose a © Wi — a © u>2. Then — 102) = ^Pprime for 

some integer k. Re-writing this, we have p-^ — (wi — w 2 ) = k for some integer k. This implies 
that either -p^— e Z or ^-«a e z. Since p-^- ^ Z, e Z. Thus wj - w 2 = since 

^ p r i ru e P r & 171 e P r & 171 e P f i Tn 6 

-(-Pprime — 1) < — < -Pprime ~ L Hence W\ = W 2 - ■ 

Lemma 3: If u> ai and w b i are uniformly distributed over Fp ime , then (Sj aj0t ai + 53iA0tbi) 
mod 0A is uniformly distributed over {8A C n V(0A)}. 

Proo/: As in the proof of Lemma[T] (£\ Qi0t ai + £) { fiMbi) mod 0A = 0(Q]j «it i + Si A*bi) 
mod A) = 00(0; ctiWai © 0j PiWbi). Since u> a i and it^ are uniformly distributed over Fp prime , aiW a i 
and /SitOfej are uniformly distributed over Fp pr4me by Lemma[2] Then, i a.iW a i ©0^ ftwj,; is uniformly 
distributed over Fp ime , and 0(0, atw a i © i fiiWbi) is uniformly distributed over {A c n V(A)}, and 
finally 0<H0i a t w a i © 0j ftwbi) is uniformly distributed over {0A C D V(0A)}. ■ 

III. Channel Model 

The Gaussian Two-way Two-relay Channel describes a wireless communication scenario where two 
source nodes (Node 1 and 4) simultaneously communicate with each other through multiple full-duplex 
relays (Node 2 and 3) and multiple hops as shown in Figure [T] Each node can only communicate with 
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G> KI> KI> KD 

Yi = X 2 + Zj Y 2 = Xi + X 3 + Z 2 Y 3 = X 2 + X 4 + Z 3 Y 4 = X 3 + Z 4 

Fig. 1. The Gaussian Two-way Two-relay Channel Model. 

its neighboring nodes. The channel model may be expressed as (all bold symbols are n dimensional) 

Yj = X2 + Zi 
Y2 = Xi + X3 + Z2 
Y3 = X 2 + X 4 + Z 3 
Y 4 = X 3 + Z 4 

where Zi (i E {1, 2, 3, 4}) is an i.i.d. Gaussian noise vector with variance Nf. Zj ~ Af(0, Nil), and the 
input X, is subject to the transmit power constraint Pf. ^E(X.f Xj) < Pj. Note that since we can always 
subtract the signal transmitted by the node itself, they are omitted in the channel model expression. Also 
note the arbitrary power constraints and noise variances but unit channel gains. 

An (2 nRa , 2 nRb , n) code for the Gaussian Two-way Two-relay channel consists of the two sets of 
messages w a ,Wb uniformly distributed over M a := {1,2, ••• ^2 nRa } and Mb := {1,2,- •• ,2 nRh } 
respectively, and two encoding functions X™ : M a — > K™ (shortened to X 4 ) and X% : Mb — > R n 
(shortened to X4), satisfying the power constraints Pi and P4 respectively, two sets of relay functions 
{fk,j}J—i (k — 2, 3) such that the relay channel input at time j is a function of the previously received 
relay channel outputs from channel uses 1 to j — 1, Xk.j = fk.j(Yk.i, ■ ■ ■ >5"fe,j-i)> an d finally two 
decoding functions gi : y[ l x M. a — > Mb and g A : y% x Mb — > M a which yield the message estimates 
Wb '■= 9x(¥i, w a) and w a := g2(Y^,Wb) respectively. We define the average probability of error of the 
code to be P n , e := 2n(R l +Rb) J2w a eM a ,w b eM b Pr{(w a ,w fe ) ^ (w a ,w b )\(w a ,w b ) sent}. The rate pair 
(R a ,Rb) is then said to be achievable by the two-relay channel if, for any e > and for sufficiently 
large n, there exists an (2 nRa , 2 nRb , n) code such that P n>e < e. The capacity region of the Gaussian 
Two-way Two-relay channel is the supremum of the set of achievable rate pairs. 

IV. Lattice Codes in the BC Phase of the Two-way Relay Channel 

The work (3), |4) introduces a two-phase lattice scheme for the Gaussian Two-way Relay Channel 
1 H 2 f> 3, where two user nodes 1 and 3 exchange information through a single relay node 2 (all 
definitions are analogous to the previous section): the Multiple-access Channel (MAC) phase and the 
Broadcast Channel (BC) phase. In the MAC phase, the relay receives a noisy version of the sum of 
two signals from both users as in a multiple access channel (MAC). If the codewords are from nested 
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lattice codebooks, the relay may decode the sum of the two codewords directly without decoding them 
individually. This is sufficient for this channel, as then, in the BC phase, the relay may broadcast the 
sum of the codewords to both users who may determine the other message using knowledge of their 
own transmitted message. In the scheme of [3|, the relay re-encodes the decoded sum of the codewords 
into a codeword from an i.i.d. random codebook while Q uses a lattice codebook in the downlink. 

In extending the schemes of |3), Q to multiple relays we would want to use lattice codebooks in 
the BC phase, as in |4|. This would, for example, allow the signal sent by Node 2 to be aligned with 
Node 4's transmitted signal (aligned is used to mean that the two codebooks are nested) in the Two-way 
Two-relay Channel: lf)2f>3f)4 and hence enable the decoding of the sum of codewords again at 
Node 3. However, the scheme of pj is only applicable to channels in which the SNR from the users to 
the relay are symmetric, i.e. — Sl _ \ n this case the relay can simply broadcast the decoded (and 
possibly scaled) sum of codewords sum without re-encoding it. Thus, before tackling the Two-way Two- 
relay channel, we first devise a lattice-coding scheme for the BC phase in the Two-way Relay Channel 
with arbitrary uplink SNRs ^ j^- . 



A. Lattice codes for the BC phase for the general Gaussian Two-way Relay Channel 

In this section, we design a lattice coding scheme for the BC phase of the Gaussian Two-way (single- 
relay) Channel with arbitrary uplink SNR - i.e. not restricted to symmetric SNRs as in H). A similar 
technique will be utilized in the lattice coding scheme for Two-way Two-relay Channel in Section [V] 

The channel model is the same as in (3): two users Node 1 and 3 communicate with each other 
through the relay Node 2. The channel model is expressed as 

Yj = X 2 + Z 1 

Y 2 = Xi + X3 + Z 2 

Y3 = x 2 + z 3 

where Z\ (i g {1,2,3}) is an i.i.d. Gaussian noise vector with variance Nf. Zi ~ A/"(0, A^I), and the 
input Xi is subject to the transmitting power constraint Pc ^E(X.f'Xi) < Pi. Similar definitions of 



codes and achievability as in Section III are assumed. 

We will devise an achievability scheme which uses lattice codes in both the MAC phase and BC 
phase. For simplicity, to demonstrate the central idea of a lattice-based BC phase which is going to be 
used in the Two-way Two-relay Channel, we do not use dithers nor MMSE scaling as in |3j-|5[^] 

'Dithers and MMSE scaling allows one to go from achieving rates proportional to log(SNR) to log(l + SNR). However, 
we initially forgo the "1+" term for simplicity and so as not to clutter the main idea with additional dithers and MMSE scaling. 
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We assume that Pi = N 2 p 2 and P3 = p 2 where p £ K and N E Z. This assumption will be 
generalized to arbitrary power constraints in the next section. We focus on the symmetric rate for the 
Two-way Two-relay Channel, i.e. when the coding rates of the two messages are identical. 

Codebook generation: Consider the messages w a , Wf, <E Fp prime = {0, 1,2,..., Pprime — !}■ Pprime 
is a large prime number such that P pr ime = [2 nRsym ], where R sym is the symmetric coding rate and [ 
denotes rounding to the nearest prime (P pr ime — [2 nRsym ] — > 00 as n — > 00 since there are infinitely 
many primes). The two users Node 1 and 2 send the codewords X x = Npt a = Npcf>(w a ) and X 2 = 



pt b = p<p{wb) where <j>{-) is defined in ([3]) in Section II-B with the nested lattices A C A c . Notice 
that their codebooks are scaled versions of the codebook C\ c y. The symmetric coding rate is then 

R —A loe ^ L 

In the MAC phase, the relay receives Y2 = Xi + X3 + Z2 and decodes (Npt a + pt b ) mod NpA 
with arbitrarily low probability of error as n — > 00 with rate constraints 



log 



2 l0g lt 



i + 



according to Lemma |4] 

Lemma 4: Let X a = a0t a = a9(f)(w a ) € {a6A c n V(a9A)} and X b = 9t b = 9<j){w b ) £ {9A C n 

,, a e Z+, 9 e R + and 



II-B 



with 



V(9A)} where w a ,Wb € Fp prime , a e Z+, 9 6 M + and 0(-), A C A c are defined as in Section 
R := i log ) ■ F rom the received signal Y = X a + X5 + Z where Z ~ A/"(0, cr^I) one may decode 
(a9t a + 9tb) mod a9A with arbitrary low probability of error as n — » 00 at rates 



1 (7 2 (a6lA) 
1 <7 2 {9AY 



1 + 



Proof: The proof generally follows p), except that for simplicity, we do not use dithers or MMSE 
scaling of the received signal in our scheme. The receiver processes the received signal as 

Y mod a9A = X a + X 6 + Z mod a9A 
= a9t a + 9t b + Z mod a9A 

To decode (a9t a + 9tb) mod a9A, the effective noise is given by Z with variance a 2 rather than the 
equivalent noise after MMSE scaling as in J3J. All other steps remain identical. The effective signal- 
to-noise ratios are SNR a — a < - ae 2 Aa ^ and SNR b = a resulting in the given rate constraints. 
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for the 1 <— 2 direction. While the rate constraint for the direction 2 — > 3 is as large as expected, the 
rate constraint for the direction 1 2 does not fully utilize the power at the relay, i.e. the codeword t& 
appears to use only the power P2/N rather than the full power P%. One would thus want to somehow 
transform the decoded sum (Npt a + pt b ) mod NpA such that both t Q and t b of the transformed signal 
would be uniformly distributed over V(-\ZP~2~A). Notice that the relay can only operate on (Npt a + pt b ) 
mod NpA rather than Npt a and pt b individually. 

Re-distribution Transform: To alleviate this problem we propose the following "Re-distribution 
Transform" operation which consists of three steps: 

1) multiply the decoded signal by N to obtain N((Npt a + ptb) mod NpA), 

2) then take mod A to obtain 



according to the operation rule in |2]). Notice that (N-\/P^t a + y/P^tb) mod y/P^A is uniformly 



dimensional lattice in order to gain some intuition and insight (though we operate in n-dimensions in our 
achievability proof). In this simple example, Node 1 sends 2pt a E {2pA c n V(2pA)} = {—2, — 1, 0, 1, 2} 
and Node 3 sends pt b e {pA c n V(pA)} = {-1, -1/2, 0, 1/2, 1}. The relay decodes 

(2pt a + pt b ) mod 2pA e { P A C n V(2pA)} = {-2, -3/2, -1, -1/2, 0, 1/2, 1, 3/2, 2, 5/2} 

and transforms it into (Apt a + 2pt b ) mod 2pA <E {2pA c n V(2pA)} = {-2, -1, 0, 1, 2} according to the 
three steps in the Re-distribution Transform. Notice |{2pA c n V(2pA)}| < \{pA c n V(2pA)}|, and hence 
the transformed signal is thus easier to forward. 



The relay broadcasts X 2 = {Ny/p2t a + \fP~2t b ) mod y/P^A. Notice that (NyfP~2t a + y/P^t b ) 
mod i/P^A is uniformly distributed over {^/P^Ac n V(\/ftA)}, and so its coding rate is R sym . Node 



N((Npt a +pt b ) mod NpA) mod A = (iV 2 pt a + Npt b ) mod A^A 



according to the operation rule in Q, and finally 
3) re-scale the signal to be of second moment P2 as 




Np 



A 



distributed over {y/PiA c n V( v / i^A)} by Lemma |] 
The three steps of the Re-distribution Transform procedure are illustrated in Figure |2]for a simple one 
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User 1 sends: 2pt a S {2pA c n V(2pA)} 



-2-10 1 2 

User 2 sends: ph e {pA c n V(pA)} 
• • • * • 

-1 -1/2 1/2 1 

Relay decodes: {2pt a + pt b ) mod 2pA e {pK n V(2pA)} 



-2 -3/2 -1 -1/2 1/2 1 3/2 2 5/2 
Re-distribution Transform: 2({2pt a + pt b ) mod 2pA) = (4jjt a + 2ptt) mod 4pA 



-4 -3 -2 -1 1 2 3 4 5 

(Apt a + 2pt b ) mod 4pA mod 2pA = (Apt a + 2pt b ) mod 2pk 



-2-10 12 

\/Pr / / / 

Scaling: ^— -((4pt„ + 2pt b ) mod 2pA) = (2y/P a t a + V p fi*i) mod VPrA 



Fig. 2. Re-distribution Transform illustration for a one-dimensional lattice. 



1 and Node 3 receive Yi = X2 + Zi and Y3 = X2 + Z3 respectively and, according to Lemma [5] may 
decode (N^/P^t a + \/t\t b ) mod ^/P^A at rate 



1, Pi 

- log 

2 8 Ni 

1, P 2 

- log — 

2 B JV 3 



Lemma 5: Let X = f?t 
are defined as in Section 



II-B 



<j>{w) e {0A C n V(0A)} where 6 e M+, w e F P;pr . me and <j>\ 
with iZ = i log r/(X7 ■ From the received signal Y = X 



),AC A c 
Z where 



Z ~ A/"(0, erf I), one may decode Ot with arbitrary low probability of error as n — > 00 at rate 

i?< 2 log^ r . 

Proof: The proof generally follows J5] Theorem 5], except that we do not use dithers nor MMSE 
scaling in our scheme. The receiver processes the received signal as 



Y mod A = X + Z mod A 
= t + Z mod A 
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To decode t, the effective noise is Z with variance a 2 rather than an equivalent noise after MMSE as in 
jsj Theorem 5]. All other steps are identical. The effective signal-to-noise ratio is thus SNR = a f£f^ , 
and we obtain the rate constraints as in the lemma statement. ■ 
Nodes 1 and 2 then map the decoded (iVPfjt a + P^t/,) mod PrA to Nw a Wb by Lemma [T] 
With side information w a , Node 1 may then determine Wb', likewise with side information Wb, Node 2 
can obtain Nw a and then determine w a by Lemma [2] 



V. Two-way Two-relay Channel 

We first consider the full-duplex Two-way Two-relay channel where every node transmits and receives 
at the same time. For this channel model we first obtain an achievable rate region for special relationships 
between the power constraints at the four nodes in Theorem |6| and the related Lemma [7] We then use 
these results to obtain an achievable rate region for the general (arbitrary powers) Gaussian Two-way 
Two-relay channel, which we show is to within i log(3) bits/s/Hz per user of the symmetric rate capacity 
in Theorem [8] 

Theorem 6: For the channel model described in Section III if Pi — p 2 , P2 = M 2 q 2 , P3 = N 2 p 2 
and P4 = q 2 , where p,q £ K + and M, N £ Z + the following rate region 

~1 



R a ,Rb < min 



2 & \N 2 



^ iog (S 



2 & \N 3 



1+ r 



1, fP 3 



1 + r 



log 



log 



P 3 

P2 



1 + N 



(4) 



is achievable using lattice codes. 
Proof: 

Codebook generation: We consider the good nested lattice pair A C A c with corresponding codebook 
Ca c ,v = {A c n V(A)}, and two messages w a ,w b £ ^p prime = {0, 1, 2, . . . , P prime - 1} in which 
Pprime is a large prime number such that P pr i me = [2 nR3ym ] {R sym is the coding rate). The codewords 
associated with the messages w a and Wf, are t a = <f>(w a ) and = <fi(wb), where the mapping <fi(-) from 
IF Ppri mc to ^A c ,v € K" defined in Q in Section 



II-B 



Encoding and decoding steps: We use a Block Markov Encoding/Decoding scheme where Node 1 and 
4 transmit a new message w a i and ivm, respectively, at the beginning of block i. To satisfy the transmit 
power constraints, Node 1 and 4 send the scaled codewords Xij = pt al = p<j>(w a i) £ {pA c n V(pA)} 
and = qtbi — q<j>(wbi) G {qA c D V(qA)} respectively in block i. Node 2 and 3 send 'X.2% and 
X3j, and Node j (j = {1,2,3,4}) receives Yji in block i. The procedure of the first few blocks (the 
initialization steps) are described and then a generalization is made. We note that in general the coding 
rates R a for w a and Rb for Wb may be different, as long as R aym = max(P a , Rb), since we may always 
send dummy messages to make the two coding rates equal. 
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O KD 



P 2 = M 2 g 2 



p 3 = a V 



<I> 



P^q 2 



Block 1: 

Send: pt al 
Decode: 
Block 2: 

Send: Pta2 
Decode: 

Block 3: 
Send: pt o3 

Decode: Wf>i (pt a3 + ATpt^ + A/ATpt^) mod NpA 



Block 4: 

Send: pt a4 

Decode: w b2 

Block 5: 

Send: P*a5 



Mgt ai 

(pt„ 2 +JVpt M ) modATpA 
(Mgt a2 + NMqt b i) mod MgA 



(Mgt a3 + iVM g t 62 + 
NM 2 qt al ) mod A/gA 

(pt a4 + Npt b3 + A/Afpt a2 
+MA" 2 pt 61 ) mod AfpA 

(Mgt a4 + NMqt b3 + ATM 2 gt a2 
+AT 2 M 2 gt 6 i) mod MqA 



Decode: w b3 (ptaS + Npt b4 + MNpt a3 + A/A" pt b2 

+M 2 N 2 pt al ) mod NpA 



gtbi 

Npt bl 

(qt b2 + Mqt a i) mod MqA 

{Npt b2 + MNpt a i) mod NpA 
(qt b3 + Mqt a2 + NMqt bl ) mod MqA 



(Npt b3 + MNpt a2 + 
MN 2 pt bl ) mod NpA 

{qt b4 + Mqt a3 + NMqt b2 + 
NM 2 qt al ) mod MqA 

{Npt b4 + MNpt a3 + MN 2 pt b2 
+M 2 N 2 pt al ) mod NpA 

(qt b5 + Mqt ai + NMqt b3 + NM 2 qt a2 
+N 2 M 2 qt bl ) mod MqA 



qt b 4 

w a2 

9*65 

W a3 



Block i: 

Send- ptal Wl^-V + NMqt b(l _ 2) + NM 2 qt a(t _ 3) 
+ ... + Ar(»-i)/2 A f( i - 1 )/ 2 gt 61 ) modAfgA 



Decode: w b{i _ 2) ( ptai + Npt^.^ + MNpt a(i _ 2) + ... 

+M (i-i)/2 Ar ( i -i)/2 pta ^ mod NpA 



(Apt^.,.!) + A/Apt a(i _ 2) + MN 2 pt b(l _ 3) qtu 
+ ■■■+ M^-WM'-V^pta!) mod NpA 

(qtu + Afgt o(< _i) + NMqt b(i _ 2) + ... w a{i _ 2) 
+N (i-D/2 M ( i - 1)/2 qt bl ) mod MqA 



Fig. 3. Multi-phase Block Markov achievability strategy for Theorem [6] 



Block 1: Node 1 and 4 send new codewords Xn = pt a i and X41 = qt^i to Node 2 and 3 
respectively. Node 2 and 3 can decode the transmitted codeword with vanishing probability of error if 



(5) 
(6) 



according to Lemma [5] 
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Block 2: Node 1 and 4 send their respective new codewords X12 = pt a 2 and X42 = qt b 2, while 
Node 2 and 3 broadcast X22 = Mqt a \ and X32 = Npt b \ received in the last block. Note they are scaled 
to fully utilize the transmit power P 2 = M 2 q 2 and P 3 = N 2 p 2 . Node 2 receives Y 2 2 = X12 + X32 + Z22 
and decodes (pt a 2 + Npt b i) mod NpA with arbitrarily low probability of error if R satisfies <|3j and 



log 



N 2 



Similarly, Node 3 can decode (qt b 2 + Mqt a i) mod MqA subject to |6]i and 



Rsym ^ 







+ 


\ log 













(7) 



(8) 



Block 3: 

Encoding: Node 1 and 4 send new codewords as in the previous blocks. Node 2 further processes 
its decoded codewords combination according to the three steps of the Re-distribution Transform 
from previous block as 

f2„ 



(N((pt a2 + Npt bl ) mod NpA)) mod NpA = {Npt a2 + N 2 pt bl ) 

= {Npt a2 + N 2 pt bl ) 
including scaling to fully utilize the transmit power P2 = M 2 q 2 as 



mod N 2 pA mod NpA 
mod NpA 



Mq 
~Np 



{Npt a2 + N 2 pt bl ) mod NpA = {Mqt a2 + NMqt bl ) mod MqA. 



It then broadcasts X23 = (Mqt a 2 + NMqt b i) mod MqA. Notice that since 
(Mqt a2 + NMqt bl ) mod MqA G {Af(?A c n V(MgA)} 



>-33 



according to Lemma [3] its coding rate is R sym - Similarly, Node 3 broadcasts X3 
MNpt a i) mod NpA again at coding rate R sym . 

Decoding: At the end of this block, Node 2 is able to decode (pt a 3 + Npt b 2 + MNpt a i) mod NpA 
with rate constraints |5]) and ^ according to Lemma [4] and Node 3 decodes (gt^ + Mqt a2 + 
NMqt bl ) mod MqA with constraints ^ and (^J. Node 1 decodes (Mgt a2 +A r Afgt b i) mod MqA 
sent by Node 2 as in the point-to-point channel with rate constraint 

Rsy m < ^ lo s(|r) (9) 

according to Lemma [5] From the decoded (Mqt a 2 + NMqtbi) mod MqA, it obtains w a 2® Nw b \ 
(Lemma [TJ. With its own information w a 2, Node 1 can then obtain N ® w b \ = w a 2 © Nw b \ Q w a 2, 
which may be mapped to w b i since P pr ime is a prime number (Lemma |2|. Notice P pr ime = 
\pnR sym ^ 00 as n — !• 00, so A < P prime = [2 nR ] and 
w a \ with rate constraint 





/jP2 Y 


]log( 







Z. Similarly, Node 4 can decode 

(10) 



December 7, 2012 



DRAFT 



14 



Block 4 and 5 proceed in a similar manner, as shown in Figure [3] 
Block i: To generalize, in Block i (assume i is oddj] 

• Encoding: Node 1 and 4 send new messages Xh = pt a i and X/ij = qt^ respectively. Node 2 and 
3 broadcast 

X 2l = (Afgt o(i _i) + NMqt b{l _ 2) + NM 2 qt a(l _ 3) + 7V 2 M 2 gt b(l _ 4) + • • • + N {l ~ 1)/2 M (l ~ 1)/2 qt bl ) mod MqA 
X 3l = (Nptw-i) + MNpt a(l . 2) + MN 2 P t b(l _ 3) + M 2 iV 2 pt a(l _4) + ■ ■ ■ + A/ (l ~ 1)/2 iV ( ^ 1)/2 pt al ) mod NpA. 

• Decoding: Node 1 decodes the codeword from Node 2 with rate constraint (|9]l (Lemma [5]) and maps 

it to u; (i_i) © Nw b (i_ 2 ) © NMw a (i_^ © N 2 Mw h (i-4) © • • • © A/' (<_1)/2 M (i " 1)/2_:L iy6i (Lemma 
[TJ. With its own messages ui a ; (Vi) and the messages it decoded previously {wbi,w b 2, ■ ■ ■ , w&(i-3)}> 
Node 1 can obtain N © w b (i_2) and determine Wb(i_ 2 ) accordingly (Lemma |2ji. Similarly, Node 4 



can decode w a (i_ 2 ) subject to rate constraint (10 1 



Re-distribution Transform: In this block i, Node 2 also decodes 

( P t m +iVpt b(l _ 1) +AfiVpt a(l _ 2) +7\/iV 2 ptfc(^3)+M 2 A r2 pta( l -4)+- ■ -+Af (l - 1)/2 iV (l - 1)/2 pt al ) mod NpA 

from its received signal Y 2 ^ = X^ + X 3 i + Z 2 ^ subject to rate constraints <|3j and (j7]i (Lemma HJ. 
It then uses the Re-distribution Transform to process the codeword combination as 

(N(pt ai + Npt K i-!) + MNpt a(l _ 2) + ■■■ + MWN^-Wpta! mod NpA)) mod NpA 
=(Npt ai + N 2 pt b{i _ 1} + MN 2 pt a[l _ 2) + ■■■ + M( J - 1 )/ 2 7V( l - 1 )/ 2+1 pt al mod N 2 pA) mod NpA 
=Npt ai + iVV^.i) + MN 2 pt a(l _ 2) + ■■■ + M^/ 2 N^/ 2+1 pt al mod iVpA 
and scales it to fully utilize the transmit power: 

^-(Npt ai + Ar a pt 6 ( i _ 1 ) + M/V 2 pt o(i _ 2) + • • • + M^/ 2 N^' 2+1 pt al mod iVpA) 
=Mgt Qi + NMqtf,^ + NM 2 qt a{t _ 2) + ■■■ + N^/ 2 M^' 2+1 qt al mod MqA. 

This signal will be transmitted in the next block i + Node 3 performs similar operations, decoding 

qtbi+Mqt a ^_^+NMqt b ^_ 2 )-\ h-/V( i - 1 )/ 2 M( i - 1 )/ 2 gt 6 i mod MqA subject to constraints ^ 

and (jfl}, and transforms it into Npt bi +MNpt a (i- 1 )+MN 2 pt b (i- 2 )-\ VM^-^^N^-^^ptu 

mod MqA, which is transmitted in the next block. 
Combining all rate constraints, we obtain 



R sym < min 



2 & \N 2 



^ l0g (t 



- log | — 
2 *\N 3 



l l0 Hw 2 



- log ( — 
2 & I iV 4 



^ l0g (f 



2 For i even we have analogous steps with slightly different indices as may be extrapolated from the difference between Block 
4 and 5 in Fig. |3J which result in the same rate constraints. 
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Assuming there are I blocks in total, the final achievable rate is ^j^R sy m, which, as / — > oo, approaches 
Rsym and we obtain Q. ■ 

In the above we had assumed power constraints of the form P\ = p 2 , P2 = M 2 q 2 , P3 = N 2 p 2 
and P4 = q 2 , where p, q E M. + and M,N € Z + . Analogously, we may permute some of these power 
constraints to achieve the same region as follows: 

Lemma 7: The rates achieved in Theorem [6] may also be achieved when Pi = N 2 p 2 , P3 = p 2 
and/or P 2 = q 2 , P 4 = M 2 q 2 . 

Proof: The proof follows the same lines as that of Theorem [6] and consists of the steps outlined in 
in Figure [5] in Appendix [A] In particular, since the nodes have different power constraints the scaling of 
the codewords is different. However, as in the previous Theorem, we only need Xi and X3 to be aligned 
(nested codebooks), and X2 and X4 to be aligned. As in Theorem [6] the relay nodes again decode the 
sum of codewords, perform the Re-distribution Transform, with a new power scaling to fully utilize their 
transmit power, and broadcast the re-distributed sum of codewords. ■ 

Theorem [6] and Lemma [Tjboth hold for powers for which P1/P3 and/or P2/P4 are either the squares 
of integers or the reciprocal of the squares of integers. However, these scenarios do not cover general 
power constraints with arbitrary ratios. We next present an achievable rate region for arbitrary powers, 
obtained by appropriately clipping the power of the nodes such that the new, lower powers are indeed 
either the squares or the reciprocals of the squares of integers. We then show that this clipping of the 
power at the nodes does not result in more than a | log(3) bits/s/Hz loss in the symmetric rate. 

Theorem 8: For the Two-way Two-relay Channel with arbitrary transmit power constraints, any rates 
satisfying 



R a ,Rb < R a , 



chievabLe 



■ max mm 

P'<P i , I ^-=N 2 or 1 ,^1=M 2 or -A* 

1 — p!* N 2 ' P{ M 2 



log 





+ 


jlogl 













-log — 
2 s \N 4 



1+ r 



1 



log 



1+ r 



- log — 
2 & \N 2 



log 



iVi 



(11) 



for some N, M € Z + and i € {1,2,3, 4}, are achievable. This rate region is within | log 3 bit/Hz/s per 
user from the symmetric rate capacity. 

Proof: Since P1/P3 and/or P2/P4 are in general neither the squares or the reciprocals of the squares 
of integers, we cannot directly apply Theorem [6] Instead, we first truncate the transmit powers to P/ 
such that P! satisfy either the constraints of Theorem [6] or Lemma [7] The achievable rate region then 
follows immediately for the reduced power constraints P[. For example, if P\ — 1 and P3 = 3.6, we 
may choose P[ = 0.9 and P3 = 3.6 , so that P$/P[ = 2 2 . Optimizing over the truncated or clipped 
powers yields the achievable rate region stated in Theorem [8] 

An outer bound to the symmetric capacity of the AWGN Two-way Two-relay Channel is given by 
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the minimum of the all the point-to-point links, i.e. 

N 3 



Ra, Rb < R 



min C 



ft 

N 2 



,c 



c 



ft 



,c 



Pi 

N 3 



,C 



ft 
N 2 



,C 



ft 

iVi 



To evaluate the gap between (111 and (IT2b: 



Rachievable "i" ^ ^§ ^ 



(a) 



max 

Pi<Pi,Pi<P 3 S 



=N 2 1 



loe 



: max 

Pi<Pl.P^<P 3 ,^r=N 2 or 



max mm 



EL 

N 2 



log 



W / . /I, (Wt 
> max I mm \ - log I 



(c) 

> min 



>Ro 



log 1 



ft_ 
JV 2 



1 

'2 



log 



log 1 



ft 

iV 4 



3P{ 



log 



loe 



1 



P' 

zJL 

Na 



loe 



P' 

zJL 

No 




- log 

2 8 V ^4 



1, /3P^ 



log 3 



3P| 

A 2 



log 3 



log 1 



ft 

A 2 



(12) 

(13) 
log 3 
(14) 

(15) 
(16) 

(17) 
(18) 



The first equality (a) follows from an assumption that WLOG, one of | log (j^) , \ log 
or ilog^-^0 is the tightest constraint. The first inequality (b) follows since the rates achieved 
by the optimized powers must be larger than those achieved by one particular strategy that meets the 
constraints - in this case the strategy that yields the P* which we construct in Appendix |B| Inequality (c) 



follows from the fact that 2P* > Pi, as also shown in Appendix 



Finally, we bound i log 



3P[ 

N 2 



with 



i log (l + as follows: If § > 2, it follows that g > 1, and hence \ log > \ log (l + ^) 

Otherwise, jlog^l + ^^ < |log3. Similarly, we may bound |log^^-^ and ^log^^f-) 
§ !og (l + and § log (l + j^J respectively. 



with 



Remark 1: We achieve a constant gap to capacity for the symmetric rate of | log(3) bit/s/Hz. The 
only other scheme we are aware of that has been shown to achieve a constant gap is noisy network 
coding pO) which achieves a larger gap of 1.26 bits/s/Hz to the capacity. The improvement in rates 
of the proposed scheme may be attributed to the removal of the noise at intermediate relays (it is a 
lattice based Decode and Forward scheme) without the sum-rate constraints that would be needed in 
i.i.d. random coding based Decode-and-Forward schemes. 



VI. Extensions to half-duplex channels and more than two relays 

We now extend our results to half-duplex channels and to two-way relay channels with more than 
two relays. 
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A. Half-duplex 

The proposed lattice coding scheme may be generalized to channels with half-duplex nodes, i.e. in 
which a node may either transmit or receive at a given time but not both. The scheme is illustrated in 
Figure |4] where we see that the proof generally mimics the full-duplex case (Theorem [6] Lemma [7] and 
Theorem [8), but that each phase (or block) in the full-duplex case is divided into two phases / blocks 
in the half-duplex case, as nodes may not transmit and receive at the same time. Thus, one may achieve 
half the rates as in the full duplex case, i.e. 

1 

R a ,R b < Rhalf-duplex = ^achievable, 

where R ac hievabie is expressed in Theorem [8] 

B. More than two relays 

This lattice coding scheme may also be generalized to more than two relays. For example, for the 
Two-way Three-relay Channel with five nodes: 1«2«3«4h5, we may apply the same strategy 

by aligning the codewords between 1 and 3, 3 and 5, and 2 and 4 (i.e. nest the corresponding codebooks). 

p' p' p' 

To align or nest the codebooks, one may truncate the transmit powers so that -wr and are either 
squares, or the reciprocals of squares of integers. By extending our Block Markov strategy, the final 
achievable rate region would be: 



Ra,Rb < 



max 

4 = N 2 or 1 
3 = M 2 01 -i, 
4 = K 2 or A, 



mm 

k = {1, 2, 3, 4}, 
j = {2,3,4,5} 



log 



N k+ _ 



log 



P' 

Ni-i 



where M, N,K G Z+ and i = {1, 2, 3, 4, 5}. 



VII. Conclusion 

We proposed a lattice coding scheme for the AWGN Two-way Two-relay Channel (1 o 2 f> 3 f> 4) 
which achieves within | log 3 bit/Hz/s from the symmetric rate capacity. This scheme may be generalized 
to half-duplex nodes, and two-way channels with more than two relays. In our scheme, each relay decodes 
a sum of codewords as all transmitted signals are properly chosen lattice codewords, performs the "Re- 
distribution Transform" which maps the decoded lattice point to another so as to fully utilize its transmit 
power, and broadcasts this transformed, scaled, lattice codeword. All decoders are lattice decoders and 
only a single nested lattice codebook pair is needed in our scheme. Note that this work results in a 
symmetric rate region which means the rates of both directions have to satisfy the constraints of all the 
links. One interesting open question is whether one can derive an achievability scheme which would 
result in an asymmetric rate region (where rates of each direction are only constrained by links for that 
direction and the two directions would thus not interfere), and whether this would result in better rates. 
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© KD ®- © 

pt a i pt al qt b i qt bl 

Q -©- Q Q 

pt a2 (pt„2 + Npt bl ) mod NpA Npt bl 

© © KD- © 

Mqt al (qt b2 +Mqt al ) mod MqA qt b2 

©- Q -©- Q 

w bl (Mqt a2 +NMqt bl ) mod MqA (qt b3 + Mqt a2 + qtt/l 

NMqt bl ) mod MqA 

© M©- © H® 

pt a i (Pt«3 + ^,,2+ (Npt b2 + MNpt al ) mod NpA w al 

MNpt al ) mod NpA 

G> © KI> © 

w b2 (Mqt a3 + NMqt b2 + (qt u + Mqt a3 + NMqt b2 + qtbi 

NM 2 qt al ) mod MqA NM 2 qt al ) mod MqA 

Q -©- © -® 

pt al (ATptj^.i) + MA r pt a( ,_2 ) + MJVV b( ,_ 3) 

+ ■ ■ ■ + M^-WNt'-Wpta!) mod JVpA 

(pt a , + Afjrtwj-H + MW Pta(,-2) + ■ ■ ■ 

w a(i-2) 

+ M <i - 1 ' /2 7V ( - 1)/2 ptai) mod JVpA 
©- © Q 



(Afgt a(i _ 1) + NMqt b{l _ 2) + NM 2 qt a( ,_ X) 
+ ■■■+ A r< '- 1)/2 M < '- 1 ' /2 9 t 61 ) mod MqA 

(qtu + Mqtad-t) + NMqt Hi - 2) + ... 
+A r( ! -l)/2 M (,-l)/2 (3tw j mod MgA 



Fig. 4. Half-duplex case. 
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P 1 = N 2 p 2 



P 2 = M 2 q 2 



<2> 



Block 1: 

Send: Npt al 
Decode: 
Block 2: 
Send: Npt a2 
Decode: 
Block 3: 
Send: Npt a3 
Decode: Wbl 

Block 4: 
Send: Npt ai 

Decode: w b 2 

Block 5: 
Send: Npt a5 

Decode: w b3 



Nptai 
Mqt al 

{Npt a2 +pt bl ) mod NpA 

(NMqt a2 + Mqt bl ) mod MqA 
(Npt a3 + pt b2 + Mpt al ) mod NpA 

{NMqt a3 + Mqt b2 + 
M 2 qt al ) mod MqA 

(Npt ai +pt b3 + NMpt a2 
+Mpt bl ) mod NpA 

(NMqt aA + Mqt b3 + NM 2 qt a2 
+M 2 qt bl ) mod MqA 

(Npt a5 + pt u + NMpt a3 + Mpt b2 
+M 2 pt„i) mod NpA 



P 3 =P 2 

<3> 

(qt b2 + Mqt a i) mod MqA 



P 4 = q 2 



gt bl 



qt b2 



(pt b2 + Mptai) modpA 9*63 

(qt b3 +NMqt a2 + Mqt b i) mod MqA w ai 

(pt b3 + NMpt a2 + qt b4 
Mpt b i) mod pA 

(qt b4 + NMqt a3 + Mqt b2 + ™a 2 
M 2 qt a i) mod MqA 

(ptu + NMpt a3 + Mpt b2 qt b5 
+M 2 pt a i) modpA 

(gt 66 + NM qt a4 + Mqt b3 + NM 2 qt a2 w a3 
+M 2 qt b i) mod MqA 



Fig. 5. Multi-phase Block Markov achievability strategy for Lemma [7] 

Appendix 

A. Multi-phase Block Markov achievability strategy for Lemma [7] 

The achievable rate region for Lemma [7] is the same as that for Theorem [6} the achievability strategy 
is essentially the same, with slight variations on the re-scaling to meet the different power constraints. 
The details of who transmits and decodes what in each phase is outlined in Figure [5] 

B. A particular choice of P* and proof that 2P* > P{ (i 6 {1, 3} J in Theorem^ 

WLOG, we assume P3 > Pi. Then, m 2 < < (m + l) 2 for some integer to £ Z + . Consider 

p* 

the following strategy for choosing P* such that ^ is a non-zero integer squared or the inverse of an 
integer squared: If m 2 < Ji^. < to(to + 1), we choose P3 = to 2 Pi and P* = Pi. Then = m p ^ > 
m^n+^P! ^ \- Thus ' 2P a ^ p 3 and p * = p i- Otherwise if m(m + 1) < ^ < (to + l) 2 , we choose 
Pf = , L^ P 3 and Pi = P 3 . Then = -, f;\ 2p > T^Vp 1 > i Thus 2P* > Pi and P,* = P 3 . 
In general, this strategy ensures that 2P* > P^. 
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